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1.  INTRODUCTION 


Tho  major  contributors  to  angular  measurement  errors  in  phased 
array  radars  include  receiver  noise,  instrument  noise,  and  target  glint. 

At  long  range,  the  errors  tend  to  be  dominated  by  noise  at  the  receiver. 

At  shorter  range  where  the  SNR  is  h’gh,  the  instr.u  lent  errors  can  pre¬ 
dominate.  At  very  short  range,  angular  error  due  to  target  glint  can  provide 
the  dominant  source  of  error.  Propagation  and  multipath  can  also  cause 
substantial  errors,  but  these  will  not  be  considered  in  this  paper. 

In  situations  where  a  radar  is  used  to  control  the  intercept  of  a  tar  got 
by  an  interceptor,  the  SNR  at  the  intercept  range  is  often  quite  high  and 
object  tracking  errors  tend  to  be  dominated  by  radar  instrument  errors. 

If,  furthermore,  the  same  radar  is  used  to  track  both  target  and  interceptor, 
fixed  or  slowly  varying  instrument  bias  errors  do  not,  to  a  first  order, 
affect  the  error  in  measuring  the  relative  position  of  target  and  interceptor. 
For  such  problems  it  is  the  randomly  varying  (i.  e.  ,  with  angle)  component 
of  the  angular  instrument  errors  which  is  of  greatest  concern. 

Attention  here  is  thus  confined  to  the  randomly  varying  component  of 
phased  array  angular  instrument  errors.  The  effect  of  these  errors  on 
tracking  accuracy  depends  both  on  their  magnitude  and  spatial  correlation 
properties,  i.  e.  ,  the  manner  in  which  these  errors  change  with  beam 
pointing  angle.  In  many  analyses  these  instrument  errors  are  assumed  to 
be  statistically  uncorrclated  from  pulse  to  pulse.  In  practice  these  errors 
will  tend  to  decorrelate  at  a  rate  which  depends  on  the  angular  rate  of  the 
object  being  tracked.  The  assumption  that  these  errors  arc  statistically 
unco rrelateii  can  lead  to  a  serious  underestimation  of  angular  tracking 
errors  (and  hence  interceptor  miss  distance). 

The  principal  contributors  to  the  random  angular  instrument  errors 
considered  in  the  analysis  arc:, 

a.  Finite  computation  accuracy  in  the  beam  steering  computer, 
including  word  truncation  in  row  and  column  adders 
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b.  Truncation  of  steering  commands  in  the  digitally  controlled 
array  phase  shifters 

c.  Phased  array  element  differential  phase  errors. 

Those  phase  errors  in  the  antenna  aperture  illumination  produce  errors  in 
the  position  of  the  null  in  the  monopulsc  difference  pattern.  An  error  can 
also  arise  from  an  error  in  the  slope  of  the  mono  pulse  difference  pattern, 
but  it  will  be  assumed  here  that  the  target  being  tracked  is  very  close  to 
the  null  of  the  difference  pattern  so  that  this  error  can  be  safely  neglected. 

The  manner  in  which  phase  errors  influence  angular  measurement 
errors  depends  on  the  statistical  properties  of  these  phase  errors  across 
the  array  face,  A  spatially  periodic  behavior  in  the  errors  can  produce 
very  large  angular  measurement  errors  at  certain  pointing  angles,  (See, 
for  example,  Sec,  7.2  of  [l], )  It  is  well  known  [l],  [2]  that  the  use  of 
space  feed  collimation  or  the  addition  of  a  random  phase  shift  (appropriately 
compensated  in  the  beam  steering  command)  to  each  element  will  tend  to 
smear  out  these  spatial  periodicities  so  that  the  angular  measurement  errors 
do  not  become  very  large  in  certain  directions. 

Additionally,  certain  techniques  can  be  employed  in  the  radar  to  ensure 
that  these  errors  become  largely  decorrelated  from  pulse  to  pulse.  These 
include  the  use  of  angle  dithering,  frequency  dithering  and  phase  dithering. 
The  effectiveness  of  these  dithering  techniques  has  been  analyzed  with  the 
aid  of  a  specially  developed  computer  program  for  modeling  phased  array 
angular  instrument  errors.  The  results  are  illustrated  for  a  typical  phased 
array  configuration. 
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11.  METHOD  OF  ANALYSIS 


Beam  steering  in  a  planar  phased  array  antenna  is  achieved  by 
establishing  a  linear  phase  gradient  across  the  aperture  which,  in  turn, 
provides  the  desired  direction  of  wave-front  propagation.  As  is  expected, 
amplitude  and  phase  errors  across  the  aperture  face  will  cause  errors  in 
both  the  transmitted  wave  front  direction  and  the  measured  angle-of- arrival 
of  reflected  energy.  Angle-of-arrival  of  energy  reflected  from  an  object 
is  determined  by  processing  both  the  output  of  the  sum  beam  and  the  mono¬ 
pulse  difference  output  [3].  It  has  been  shown  by  Nester  [4]  that  when 
amplitude  and  phase  errors  across  the  aperture  are  small,  angle  measure¬ 
ment  errors  are  given  to  a  first  order  by  a  weighted  sum  of  the  phase  errors 
in  the  monopulse  difference  pattern  illumination  function.  In  particular, 
Nester's  analysis  shows  that  angle  measurement  errors  are  independent,  to 
a  first  order,  of  sum  pattern  amplitude  and  phase  errors  and  difference 
pattern  amplitude  errors. 

The  phase  error  at  each  array  element  in  the  monopulse  difference 
pattern  illumination  function  consists  of  a  bias  component,  i.e.  ,  insertion 
phase  error,  and  a  variable  error  which  depends  on  the  commanded  phase. 
The  insertion  phase  error,  as  Nester  shows,  gives  rise  to  a  fixed  instrum.ent 
bias  error.  Since  instrument  bias  errors  are  of  no  concern  here,  insertion 
phase  errors  will  be  ignored  in  this  paper.  The  variable  phase  errors  are 
due  to  finite  computation  accuracy  in  the  beam  steering  computer,  truncation 
of  steering  commands  in  the  digitally  controlled  phase  shifters,  and  element 
differential  phase  errors.  The  latter  errors  are  associated  with  the  bit 
setting  of  each  phase  shifter  and  will  be  defined,  along  with  the  other  error 
sources  in  Sec.  II-C. 

A,  Phased  Array  Geometry 

The  antenna  face  is  assumed  to  consist  of  a  regular  planar  lattice  or 
elements.  A  typical  array  geometry  is  shown  in  Fig.  1.  As  shown,  x  and  > 
are  orthogonal  car  esian  coordinates  across  the  aperture  face,  and  the  /. 


coordinate,  normal  to  the  face  aperture,  is  chosen  to  form  a  right-handed 
coordinate  system,.  The  element  spacing  is  determined  by  lattice  vectors  h 
and  V  which  are  not  necessarily  orthogonal  and  do  not  necessarily  line  up 
with  the  X  or  y  axes. 

The  beam  dirfe..Hon  is  determined  by  the  angles  o  and  p,  where  a  is  the 
angle  between  the  beam  and  the  ya  plane,  and  p  is  the  angle  between  the  beam 
and  the  xz  plane,  as  shown  in  Fig.  1.  It  is  most  convenient  to  express  point¬ 
ing  angle  by  sine  coordinates  u  =  sino  and  v  =  sinP.  The  array  in  Fig.  1  is 
shown  as  circular  with  an  equilateral  triangular  grid  structure,  but  the  methofi 
of  analysis  is  not  restricted  to  this  case.  The  geometry  of  Fig.  1  is  assumed 
for  all  the  cases  analyzed  in  this  paper,  however. 

The  position  of  a  given  element  in  the  lattice  is  characterized  by 
indices  m  and  n,  and  the  vector  displacement  of  the  mnth  element  in  the 
aperture  face  is  given  by 


■  mn 


m]i  +  n_y 


(1) 


The  X  and  y  components  of  the  vector  are  simply 


■‘mn  =  "’'‘x  +  "-x 


y  -  mp.,  +  ni' 

'  mn  '^y  y 


1  i  I 


where  a  and  u  are  the  x  and  y  components  of  u  and  v  and  v  are  the  x  and 
X  y  X  y 

y  components  of_v. 

In  the  lattice  structure,  the  array  elements  can  also  be  put  on  a 
rectangular  grid  which  is  lined  up  with  the  x  and  y  axes,  as  shown  in  Fig.  1- 
The  indices  i  and  j  denote  the  column  number  and  row  number  of  the  rectan¬ 
gular  grid.  When  the  array  lattice  is  non- rectangular,  only  certain  i  j  combin¬ 
ations  are  occupied  by  phased  array  elements. 
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Error  Formulas 


B. 

The  angular  measurement  errors  are  computed  from  a  weighted  sum, 
extending  over  the  aperture,  of  phase  errors  at  each  element,  as  derived  by 
Nester  (see  Eq.  (60)  of  Ref.  [4]).  The  two  orthogonal  angle  errors, 
expressed  in  sine- space  coordinates  u  and  v  are ' 


u 


y^6  I 

xmn  mn 


m,  n 


ky6  X 

■/^°xmn  mn 
m,  n 


(4) 


\mn^mn 


m,  n 


/n,n 


(5) 


The  quantities  appearing  in  these  expressions  are  defined  as  follows: 


t  =  error  in  u 
u 


t  =  error  in  v 
V 

6^(x,y)  --  X  monopulse  difference  pattern  amplitude  weighting 
(illumination)  function  at  coordinates  x,  y 

6^(x,y)  -  y  monopulse  difference  pattern  amplitude  weighting 
(illumination)  function  at  coordinates  x,  y 


Throughout  this  analysis  pointing  angles  will  be  expressed  in  sine- space 
coordinates.  In  particular,  (he  errors  and  arc  expressed  in 
tnillisines  (n^s). 
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6  =  6  (x  ,  y  ) 

xmn  X  mn  ^mn 


6  ^  6  (x  ,  V  ) 

ymn  y  mn  ■  mn 


|(x,y)  -  phase  eri'or  at  the  co-ordinates  x,y 

I  =  |(x  ,y  ) 
mn  mn  'mn 


Ztt 

k  -  propagation  constant  =  . 


Note  that  6  is  odd  in  x  but  even  in  y,  while  6  is  even  in  x  and  odd  in  y.  It 
X  '  y  ' 

is  readily  apparent  from  (4)  and  (5;  tha'i  a  constant  phase  error  across  the 
array  face  does  not  contribute  to  the  angular  errors.  It  is  evident  that  once 
the  aperture,  i.e.  ,  shape,  lattice  structure,  etc.  ,  and  the  difference  pattern 
weighting  function  have  been  defined,  the  only  quantities  that  vary  with  scan 
angle  are  the  particular,  the  denominators  in  the  above  two  expres¬ 

sions  are  merely  normalization  constants.  It  is  of  interest,  therefore,  to 
examine  the  sources  of  the  errors  which  determine  the  I 

mn 

C .  Scan  Dependent  Phase  Errors 

Each  digitally  controlled  phase  shifter  controls  the  phase  at  each 
element  to  a  binary  submultiple  of  radians,  i.e.,  180°,  90°,  45°,  etc. 

In  practice,  each  bit  of  the  phase  shifter  introduces  a  phase  shift,  say  , 
in  the  "off"  state  and  0^^  +  2ir/2^  +  c.  in  the  "on"  state  where  i  denotes  the 
bit  number,  with  bit  t  controlling  180°,  bit  2  controlling  90°  and  so  on. 

With  all  phase  shifter  bits  in  the  "off"  position,  the  total  phase  shift  is 

This  residual  phase  is  called  the  insertion  phase.  This  insertion 
phase,  which  varies  from  element  to  element  due  to  manufacturing  tolerances, 
is  independent  of  phase  shifter  setting  and  hence  does  not  change  with  beam 
pointing  angle.  It  is  readily  seen  from  Eqs.  (4)  and  (5)  that  these  insertion 
phases  introduce  fixed  bias  errors  in  the  angular  measurements.  Since  the 
constant  insertion  phase  error  generates  only  bias  errors,  the  are  taken 
to  be  zero.  The  error  « .  represents  the  phase  error  corresponding  to  the 
ith  bit  if  it  is  "on."  These  errors  are  the  differential  phase  errors  and  result 


7 


frotn  manufacturing  tolerances.  The  statistical  properties  of  these 

differential  phase  errors  depend,  of  course,  on  the  manufacturing  pr''ccdures 

employed.  Since  these  errors  are  scan  dependent,  their  contribution  to  the 

measurement  errors,  c  and  €  ,  must  be  taken  into  account. 

u  V 

A  major  source  of  phase  error  at  each  array  elen^ent  is  due  to  truncation 
of  the  binary  phase  shift  commands  in  the  beam  steering  computer  and  at  the 
element  phase  shifters.  In  order  to  determine  how  these  truncation  errors 
affect  angle  measurement  errors,  consider  how  the  phase  shift  commands 
are  generated  for  each  array  element.  The  beam  steering  command  for  the 
ijth  element,  is  given  by  the  sum 

ij  Cl  rj  sij 


where 

=  commanded  phase  of  ith  column  for  steering  in  u 

=  commanded  phase  of  jth  row  for  steering  in  v 

=  space  feed  collimation  phase  at  coordinates  x^,  yj . 

=  0  if  no  space  feed  collimination  is  present.  The  linear  phase  gradient 
commands  for  row  and  column  steering  are  given  in  terms  of  the  angle 
coordinates  u  and  v  by 


^  .  =  kx.u  (6) 

Cl  1 

^  .  =  ky.v  (7) 


where 

X.  -  horizontal  aperture  coordinate  of  the  ith  column 
Vj  =  vertical  aperture  coordinate  of  the  jth  row. 

The  space  feed  collimation  phase  0  ..  for  the  ijth  element  is  determined  from 


where  F  is  the  focal  length  of  the  space  feed  structure. 

In  the  beam  steering  computer  the  quantities  0^.,  and  0^^^  are 

computed  to  s  bits  and  truncated  to  r  bits.  These  quantities  are  then  added 
and  the  result  truncated  to  q  bits,  which  is  the  number  of  bits  required  to 
control  the  phase  shifter. 

Figure  2  illustrates  how  these  beam  steering  operations  are  simulated 

in  the  computer  error  analysis  program  in  order  to  obtain  the  element  phase 

shift  errors  I  which  are  needed  to  determine  angular  errors  €  and  c  . 

For  any  chosen  lattice  points  m  and  n,  the  row  and  column  indices  i  and  j  are 

computed  with  the  aid  of  (2)  and  (3).  0  .  and  0  .  are  then  computed  from  (6) 

Cl  r  J 

and  (7)  using  the  commanded  pointing  angles  u  and  v  and  the  cartesian 
coordinates  of  the  elements  x.  and  y^.  If  present,  the  collimation  phase  is 
readily  calculated  from  (8).  The  row  and  column  addition  indicated  is 
performed  as  described  above  and  the  result  is  truncated  to  q  bits.  The  oox 
labeled  "compute  differential  phase  error,  etc."  represents  a  bit-testing 
function  that  senses  the  phase  shifter  bit  setting  and  determines  which  bits 
are  "on."  In  the  computer  simulation,  the  differential  phase  error  for  each 
bit  in  the  array  is  separately  generated  with  the  aid  of  a  random  number 
generator  and  stored  for  use  by  the  program.  For  each  bit  that  is  "on,"  the 
corresponding  differential  phase  error  is  recalled  from  storage  and  added 
to  the  phase  developed  at  the  phase  shifters.  Comparing  the  phase  thus  actually 
developed  with  the  desired  phase,  results  in  the  elemental  phase  error, 
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CHOOJE  LATTICE  POINT 


PHASE  ERROR  AT  POSITION  mn 

rnn 


Fig.  2  Simulation  of  Beam  Steering  Computer  Logic 
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Reference  to  P'ig.  2  shows  that  the  collimation  angle  0  ..  has  to  be 
computed  or  recalled  many  times,  typically  once  for  each  phase  shifter 
element.  It  has  been  shown  by  Hatcher  [  5]  that  this  computation  can  be 
simplified  considerably  by  expanding  this  collimation  term  as  a  sum  of  two 
terms  where  the  first  depends  only  on  x.  and  the  second  only  on  y..  This 
procedure  reduces  the  number  of  collimation  computations  to  the  sum  of  the 
number  of  rows  and  columns.  The  expression  for  the  collimation  phase  can 
be  expanded  in  a  Taylor  series;. 


0  . .  =  7-=?  (x.  +  y.  I  +  higher  order  terms  of 
’^sij  XF  \  1  ®  ,  , 

•'  '  '  even  powers  of  x.  and  y. 


2  2 

Thus  0  . .  =!  0  .  0  where  0  .  =  irx  /XF  and  0  .  -  ■nyT/K'F.  The  summation 

’^Slj  ’^si  sj  ’^si  1 

of  the  higher  order  terms,  which  are  even  powero  '  *  x^  and  weighted  by  the 
odd  difference  pattern  illumination  functions,  vanishes  identically  in  the  calcu¬ 
lation  of  and  as  defined  by  Eqs.  (4)  and  (5).  Thus  the  omission  of  these 
higher  order  terms  does  not  contribute  to  the  angular  errors  but  does  produce 
some  degradation  in  main  beam  gain  and  sidelobe  levels,  as  discussed  by 
Hatcher  [  5^  Nevertheless,  this  degradation  does  not  prove  to  be  serious  in 
most  situations  of  interest.  A  schematic  diagram  of  this  m.ode  of  operation  is 
shown  in  Fig.  3.  The  quantities  0^.,  0^^,  0^.,  and  0^^  are  computed  in  the  beatu 
steering  computer  to  an  accuracy  of  s  bits.  The  row  and  column  commands 
are  then  added  and  truncated  to  r  bits.  From  this  point  on,  the  sequence  of 
operations  in  Fig.  3  is  identical  to  that  of  Fig.  2. 

The  approximate  method  for  calculating  the  collimation  phase  as  described 
above  was  used  for  all  applicable  cases  analyzed  in  this  paper. 
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D. 


Randomized  Quantization  Levels 


A  q  bit  phase  shifter  has  uniformly  spaced  quantization  levels  (ignoring 
differential  phase  errors)  with  spacing 


The  generation  of  phase  shifter  quantization  errors  is  illustrat^'i  in  Fig  4. 

The  solid  line  represents  the  desired  phase.  The  dashed  staircase  function 
is  the  actual  phase  due  to  quantization.  The  difference  between  these  functions 
gives  a  periodic  sawtooti-  error  with  amplitude  Q/2.  This  sawtooth  is  an  odd 
function  of  commanded  phase.  An  erroj*  of  this  type  results  wherever  a  binary 
number  representing  phase  is  truncated. 

In  the  situation  where  a  linear  phase  function  is  applied  across  the 
aperture  face,  quantization  errors  will  cause  errors  which  will  be  periodic 
across  the  face  of  the  aperture.  As  noted  earlier,  spatially  periodic  phase 
errors  produce  angle  m  asurement  errors  in  certain  directions  which  can 
be  quite  large.  It  is  well  known  that  the  space  feed  collimation  correction 
tends  to  smear  out  these  periodicities  and  thereby  substantially  reduces  the 
maximum  angular  error  [l],  [2].  A  more  direct  approach  sometimes 
employed  on  corporate  fed  arrays  is  to  deliberately  vary  the  starting  position 
(but  not  the  spacing)  of  the  quantization  levels  in  a  random  way  from  element 
to  element.  This  randomization  may  be  implemented  in  some  types  of  arrays 
by  randomizing  the  length  of  cable  (or  waveguide)  which  drives  each  element 
of  the  array,  and  compensating  for  the  differential  delay  path  by  adding  an 
appropriate  compensation  phase  number  to  the  commanded  phase  for  each 
element  at  a  point  prior  to  truncation  of  the  beam  steering  command  in  the 
q  bit  phase  shifters. 

The  dramatic  improvement  in  accuracy  which  results  with  the  use  of 
randomized  quantization  levels  is  illustrated  schematically  for  a  worst  case 
situation  in  Fig.  5.  To  simplify  the  illustration,  a  one-dimensional  case 
has  been  considered  where  the  desired  phase  is  a  linear  function  of  the 
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Fig.  4  Phase  Shifter  Quantization  Errors 
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aperture  coordinate  x.  As  illustrated  in  (a),  the  beam  steering  off  the  array 
normal  is  small,  and  the  commanded  phase  causes  only  the  phase  shifters  at 
the  very  edge  of  the  aperture  to  switch  out  of  the  zero  state.  As  seen  from 
part  (b),  the  difference  between  the  commanded  and  actual  phase  functions  is 
a  linear  phase  error  term  which  will  cause  an  angular  error  approximately 
equal  to  the  commanded  steering  angle.  Part  (c)  of  the  figure  illustrates  the 
effect  of  randomized  quantization  levels.  The  horizontal  dashes  indicate  the 
location  of  the  nearest  quantization  levels  for  each  phase  shifter.  In  all 
cases,  the  commanded  phase  is  quantized  to  the  nearest  available  quantization 
level.  The  resulting  phase  errors  at  each  phase  shifter,  shown  in  part  (d)  of 
the  figure,  are  seen  to  have  a  random  character  which  will  not  give  rise  to  a 
large  angular  error. 

E.  Evaluation  of  the  RMS  Angular  Error 

A  formula  for  the  rms  angular  error  is  developed  in  this  section  for  the 
case  where  phase  shifter  phase  errors  are  spatially  uncorrelated.  The 
assumption  of  uncorrelated  errors  is  valid  only  if  some  means  of  phase 
randomization  exists,  such  as  the  method  of  randomized  quantization  levels 
discussed  in  the  previous  section.  Space  feed  collimation  also  provides  a 
form  of  spatial  decorrelation,  although  regions  of  correlation  remain,  as  will 
be  discussed  later.  The  condition  for  uncorrelated  phase  errors  is  given  by 

^mn^m'n'  "°^4^mm'^nn' 

where  the  bar  denotes  mathematical  expectation  (in  this  case  an  ensemble 

average)  and  6..  is  the  Kronecker  delta  function.  For  the  error  sources 

®  ij  2  . 

considered,  the  mean  square  phase  error  cr ^  is  the  same  for  all  elements. 
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From  Eq.  (4),  noting  that  ^  =0  (since  insertion  phase  errors  have  been 

mn  ^  _ 

set  to  zero),  it  follows  that  =  0.  The  variance  of  u,  is  then  given  by 


X)  ®xmn 

’  m,n 


IE 

'm,n 


^xmn*mn 


(9) 


An  exactly  analogous  expression  can  be  derived  for  o- 


2 

V 


These  expressions  are  easily  evaluated  numerically.  Values  of  and 
computed  for  various  difference  pattern  weighting  functions  are  given 
in  Appendix  A.  The  estimates  of  variance  from  computer  simulation 
generated  sequences  of  and  are  in  excellent  agreement  with  values 
calculated  from  the  above  formulas  in  those  cases  where  quantization  level 
randomization  or  space  feed  collimation  is  employed. 
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m.  PERIODICITY  AND  SYMMETRY  PROPERTIES 
OF  ANGULAR  ERRORS 


In  this  section  we  shall  derive  some  of  the  more  useful  periodicity 
and  symmetry  properties  obeyed  by  the  phased  array  angle  measurement 
errors.  These  properties  aid  in  the  interpretation  of  computed  error  func¬ 
tions  and  enable  the  amount  of  computation  to  be  substantially  reduced. 

These  properties  arise  from  periodic  behavior  of  the  array  lattice,  the 
particular  symmetry  properties  of  the  lattice  structure  considered,  and  the 
symmetry  and  periodicity  properties  of  phase  errors  due  to  truncation  of 
binary  steering  commands  in  each  element  phase  shifter. 

In  addition  to  the  lattice  properties  specified  by  the  lattice  vectors  ^ 
and  v_,  we  make  use  of  the  following  relationships; 

6jj(x,y)  =  6^(x, -y)  =  -6j^(-x,y)  (11) 

6y(x,  y)  =  6y(-x,  y)  =  -6y(x,-y)  (12) 

and,  as  seen  from  Fig.  4, 


m  -  -i(-0) 

(13) 

1(0)  =  e(0  +  Q) 

(14) 

where  |  is  the  phase  shifter  quantization  error  (in  the  absence  of  quantization 
level  randomization)  and  0  is  the  commanded  phase.  The  commanded  phase  at 
the  point  (x,  y)  for  a  commanded  angle  (u,  v)  is,  ignoring  collination, 

0  =  k(xu  +  yv)  .  (15) 
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PnetliRf  pvt  Unk 


With  iheso  relationships,  Eqs.  (4)  and  (5)  can  be  put  in  the  following  form, 
where  proportionality  constants  have  been  omitted: 


i  (u,  v)  ~  S  ''  6  (x  ,  y  )  I  [k(x  u  y  v)l 
u'  '  /  ^  x'  mn’  ^mn'  '  mn  ^mn  ' 


(16) 


m,  n 


€  (u,  v)  ~  "  6  (x  ,  V  )  4  tk(x  u  4  y  v)] 
v'  '  X  .«  “y'  mn  mn'  ^  '  mn  ^mn  ' 


(17) 


m,  n 


These  relationships  are  now  in  a  convenient  form  for  deriving  the  desired 
properties  of  the  error  functions. 

Let  w_  be  a  vector  with  components  (u,  v),  represents  a  vector 
pointing  angle  in  (u,  v)  space.  Periodicities  in  (u,  v)  space  are  found  by 
determining  values  of  Aw  which  leave  and  unchanged.  U  can  be 
verified  by^  direct  substitution  in  the  above  equations  that  displacements 
Aw  =  and  Aw  =  V  leave  and  unchanged,  where 


U  ^  Q(ii  x  jii  Xj;)) 

X  ii) 


k{/v^  -  (ii-v)^) 


Y  ^  Q(ti  X  X  _  Q(p^^  -  (Ji» 


Writing  €^^(u,  v)  -  it  follows  that 


€  (w  4  m’U  4  n*V  )  -  c  (w  ) 
u'—  •—  —  u  — 


(1^) 


(i'>; 


(20) 
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for  arbitrary  integers  m’  and  n*,  and  any  w .  Similarly 

€^(w  +  m'^  +  n'V^)  =  «^(w)  .  (21) 

Thus,  the  lattice  periodicity  of  the  array  leads  to  a  lattice  periodicity  in  the 
error  functions  c  and  c  characterized  by  the  lattice  vectors  U  and  V.  It 
can  be  readily  shown,  in  particular,  that  an  equilateral  triangle  lattice 
gives  rise  to  angle  error  functions  which  also  have  an  equilateral  triangle 
lattice  structure. 

Consider  now  the  symmetry  properties  of  and  in  (u,  v)  space. 
Using  (13)  together  with  (16)  readily  gives 

-v)  =  -£^(u,  v)  .  (22) 


Similarly, 


<^(-u,  -v)  =  -£^(u,  v) 


These  relations  can  be  written  more  concisely  as 


e  ( -w  )  =  -  £  (w ) 
u'  — '  u'— ' 


Equations  (24)  and  (20)  can  be  combined  to  give 


‘u(f  Ji  +  ?  i  ‘  u(r  i  -  Si ) 


(23) 


(24) 

(25) 


(26) 
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and 


vff-'  U  +  f  V  +  w)  ;  u  +  a'  V  .  w)  (27) 

lor  arbitrary  integers  m‘  and  n*.  The  error  functions  thus  have  anti¬ 
symmetry  points  about  the  basic  lattice  points  and  mid-lattice  jtoints  (i.  e.  , 
midway  between  lattice  points)  of  the  error  surface. 

Consider,  finally,  array  lattices  which  have  reflection  symmetry  about 
the  X  and  y  axes.  The  triangular  lattice  shown  in  Fig.  1  is  an  example  of 

such  a  lattice.  Then  if  we  replace  x  by  -x  and/or  y  by  -y  cvery- 

mn  ^  mn  ^mn  ’  ^mn  ’ 

where  in  the  array,  the  lattice  structural  remains  unchanged.  Replacing 
^mn  "^mn  yields,  using  (11). 


«y{vi,  v)  =  -€^(-u,  v)  .  (28) 

Replacing  y^^^  by  in  (16)  yields,  using  (11), 


v)  =  «u^u, -v)  .  (29) 

In  a  similar  manner  one  can  show 

v)  =  €^(-U,  v)  (30) 


and 


«y(u,  v)  -  -•ty(u,  -v) 


(31) 
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For  the  particular  lattice  structure  shown  in  Fig.  1,  the  vector  IJ 
is  aligned  with  the  u  axis.  For  an  arbitrary  integer  m'  it  follows  that 

t^(u  +  m'U,  v)  =  .  (32) 

Combining  this  with  (28)  gives 

.u(f'  U  +  u,v)=-.^(f'  U-u.v)  .  (33) 

One  can  show  in  a  similar  manner 

*v(r  U  v)  =  U  -  u.  v)  .  (31) 

Additional  symmetry  relationships  could  certainly  be  derived  using  similar 
methods. 

It  can  be  shown  that  pioperties  of  the  error  functions  derived  above 
still  remain  valid  if  a  space  feed  collimation  term  is  included.  These 
properties  will  be  destroyed,  however,  by  the  inclusion  of  quantization  level 
randomization  or  element  differential  phase  errors  or  both. 


IV.  ARRAY  ASSUMPTIONS  AND  FORMAT  OF  COMPUTATIONAL  RESULTS 


The  array  example  used  in  this  paper  is  a  circular  array  with  an 
equilateral  triangle  lattice  structure,  as  shown  in  Fig,  1.  The  array  parame¬ 
ters,  as  given  below,  are  representative  of  those  which  have  been  considered 
for  command  guidance  intercept  applications. 


d 

X 

(column  spacing) 

= 

0.  500  \ 

d 

y 

(row  spacing) 

= 

0. 289  X 

•^x 

(lattice  parameter) 

= 

0. 500  X 

"y 

(lattice  parameter) 

0.  289  X 

V 

X 

(lattice  parameter) 

= 

O.OOOX 

V 

y 

(lattice  parameter) 

0.578  X 

D 

(aperture  diameter) 

= 

26.4  X 

R 

(aperture  radius) 

= 

D/2  =  13.  2  X 

N 

(number  of  elements) 

= 

1891 

6 

X 

(difference  pattern 
amplitude  weighting 
function) 

sin  ttx/R 

6 

y 

(difference  pattern 
amplitude  weighting 
function) 

sin  iry/R 

F 

(focal  length  in  the 
space  feed  mode) 

= 

D 

q 

(no.  of  phase  shifter 
bits) 

- 

3 

The  boresight  beamwidth  of  this  array  is  approximately  2.6“. 
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Preceding  page  blank 


The  computer  error  analysis  program  produces  graphs  of  the  angle 

error  t  and  t  as  a  function  of  u  and  v,  and  an  autocorrelation  function  of 
u  V 

these  errors,  denoted  (Au,  Av)  and  R^(Au,  Av),  plotted  as  a  function  of 
differential  steering  angle,  Au  and  Av.  The  sample  point  spacing  for  the 
simple  planar  cuts  was  chosen  to  be  0,  25  ms.  The  computed  rms  value  of 
the  angulai  error  annotated  on  each  graph  refers  only  to  the  scan  sector 
actually  shown.  Unless  otherwise  indicated,  the  differential  phase  errors  are 
set  equal  to  zero  for  reasons  to  be  discussed  later. 

Many  of  the  plots  shcA/n  in  the  subsequent  section  are  perspective  views 
of  the  angle  errors  and  the  corresponding  two-dimensional  autocorrelation 
function.;  In  each  of  these  cases,  the  sample  point  spacing  was  chosen  to  be 
0.5  ms  in  both  u  and  v,  The  properties  of  this  perspective  plot  are  illustrated 
in  Fig,  6.  Each  of  the  perspective  plots  is  viewed  at  60 “  clockwise  from 
the  u  axis  at  a.n  elevation  angle  of  35°.  The  limits  indicated  in  the  figure 
are 


Av. 


20  ms  for  error  surfaces 


15  ms  for  autocorrelation  surfaces 


The  computation  scan  is  performed  parallel  to  the  u  axis  at  constant  v, 
similar  to  an  ordinary  raster  scan..  The  6400  point  error  surface  thus 
obtained  was  used  to  compute  the  corresponding  autocorrelation  functions. 

The  plotting  lines  are  drawn  at  the  computation  interval,  viz.,  0.  5  ms.  Edge 
effects  were  neglected  in  the  computation  of  these  autocorrelation  functions, 
which  were  generated  by  incrementing  Au  and  Av  in  one  quadrant  only.  Points 
in  the  remaining  quadrants  were  obtained  by  symmetrically  transposing  points 
about  the  origin  as  justified  by  symmetry  properties  of  the  autocorrelation 
fund  ion. 
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Fig.  6  Perspective  Views 


The  computer  error  analysis  program  is  written  in  FORTRAN  for 
use  on  the  IBM  37  0/155.  Special  truncation  and  bit  sensing  functions  that 
operate  on  the  binary  representation  of  words  are  required.  Because  of  the 
32-bit  word  structure  peculiar  to  this  machine,  reference  phases  are  com¬ 
puted  to  24  bits,  i.e.  ,  the  quantization  step  size  of  computed  phase  errors 
24 

is  2  17/2  .  The  amount  of  storage  required  for  any  particular  array  lattice 

is  given  approximately  by  3N  3  2-bit  words  plus  N(2  +  q)  16-bit  words. 
Execution  time  is  approximately  120N  psec  to  compute  both  and  for  one 
beam  position  and  3-bit  phase  shifters,  including  differential  phase  errors. 

A.  Quantization  Lobes 

In  the  absence  of  quantization  level  randomization  and  space  feed 
collimation,  the  angle  error  exhibits  large  peak  values,  called  quantization 
lobes,  near  the  array  normal  and  at  integral  multiples  of  the  error  periods 
U,  V.  These  lobes  in  the  sine  alpha  error  €  ^,  which  are  not  unlike  grating 
lobes,  are  apparent  in  the  computer  generated  plot.  Fig.  7.  Note  the  linear 
portion  of  with  slope  of  -1  in  the  region  -5  <  u  <+  5  ms.  In  this  region,  all 
phase  shifter  phase  commands  are  smaller  than  half  the  quantization  level 
spacing  Q.  Since  none  of  the  phase  shifters  switches  out  of  the  zero  state, 
the  resulting  error  is  the  negative  of  the  commanded  steering  angle.  The 
angle  at  which  the  magnitude  of  the  error  reaches  a  peak  is  approximately 
where  the  phase  shifters  at  the  edge  of  the  array  just  start  to  switch  into  the 
first  quantizing  level.  The  values  of  u  for  which  this  maximum  occurs  is 
given  approximately  by  |ul  ^  max  U  |  ^  5ms.  The  normalized  autocorrelation 
function  of  f®  shown  on  the  right  side  of  Fig.  7.  The  autocorrelation  func¬ 
tion  and  the  rms  error  in  Fig.  7  are  ‘.ot  representative  of  the  large  sector 
average  due  to  the  predominance  of  the  large  peak  in  the  sample  interval. 

Fieures  8  and  show  the  two-dimensional  error  surfaces  t  (u,v)  and  t  (u,  v), 

u  V 

fi'spect  ively .  Note  th.it  the  curv<'  in  Fig.  7  represimts  a  planar  cut  through 

the  Mirl'aci-  I  tu.vl  of  Fig.  8  normal  to  the  v  axis  and  at  v  0.  The  sym- 
u 

iiu'try  properties  as  discussed  in  Sec.  Ill  are  clearly  visible  in  both  and 
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B.  Addition  of  Space  Feed  Collimation 

The  case  presented  in  Fig.  10  is  exactly  the  same  as  the  case  presented 
in  Fig.  7  except  that  space  feed  compensation  has  been  included.  The  colli¬ 
mation  was  computed  on  a  row  and  column  basis  as  discussed  in  Sec.  II-C. 

In  comparing  this  figure  to  Fig.  7,  it  is  noted  that,  in  addition  to  a  10:1  change 
in  vertical  scale,  the  resulting  error  is  much  more  "random."  As  expected, 
the  error  function  is  still  anti -symmetric  about  the  origin,  as  discussed  in 
Sec.  III.  The  autocorrelation  function  exhibits  a  3  dB  width  of  about  1.7  ms. 
The  rms  error  is  found  to  be  0.  149  ms,  which  happen®  to  agree  well  with  the 
predicted  large  sector  rms  value  of  0.  1486  ms  calculated  from  Eq.  (9).  In 
the  case  of  the  sine  beta  error,  «  which  is  symmetric  about  the  origin,  the 
rms  error  is  0.  203  ms,  which  compares  with  a  value  of  0.  148  ms  calculated 
from  Eq.  (10).  This  apparent  discrepancy  is  due  to  the  fact  that  the  space 
feed  collimation  is  not  completely  effective  in  randomizing  phase  errors 
across  the  array  face.  Phase  errors  of  neighboring  elements  still  tend  to  be 
correlated.  Thus  the  rms  angular  errors  calculated  from  a  limited  scan 
sector  cannot  be  expected  to  conform  closely  with  errors  predicted  from 
Eqs.  (9)  and  ( 10). 

Figure  11  is  the  sine  alpha  error  surface  in  the  vicinity  of  the  array 
normal  and  shows  the  same  scanning  region  as  Fig.  8  and  9.  This  figure 
shows  that  «  ^(u,v)  is  zero  whenever  u  is  zero,  consistent  with  the  anti¬ 
symmetry  properties  discussed  in  Sec.  III.  In  fact,  careful  inspection  of 
Fig.  11  reveals  the  presence  of  other  symmetry  properties  as  discussed  in 
Sec.  III.  Figure  12  is  the  two-dimensional  autocorrelation  function  R^(Au,  Av) 
of  €  (u,  v),  i.e.,  of  the  error  surface  shown  in  Fig.  11.  The  shape  of 
R^(Au,  Av)  is  generally  similar.  This  autocorrelation  function,  as  well  as 
the  one  shown  in  Fig.  10,  shows  significant  spatial  correlation  for  differential 
angles  of  2  ms  and  beyond.  This  spatial  correlation  could  significantly  degrade 
radar  performance  for  applications  such  as  command  guidance  intercept. 
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STEERING  ANGLE,  SINE  ALPHA  (ms)  OIFF.  STEERING 

ANCLE  (ms) 


11  Perspective  View  of  the  Sine  Alpha  Error  in  the  Presence  of  Colli mation 


iTig.  12  Perspective  View  of  the  Sine  Alpha  Autocorrelation  Function 
in  the  Presence  of  Collimation 


Thus  far  the  number  of  bits  carried  in  the  beam  steering  computation 
is  large  so  that  truncation  errors  are  not  significant.  Figures  13  and  14 
illustrate  the  importance  of  row  and  column  truncation  errors.  Figure  13 
shows  the  effect  of  truncating  row  and  column  steering  commands  to  h  bits. 

The  rms  error  is  increased  noticeably  to  0.228  ms.  In  Fig,  14,  the  row  and 
column  steering  commands  have  been  truncated  to  4  bits,  and  this  results  in 
a  situation  where  the  row  and  column  truncation  errors  dominate  the  phase 
shifter  truncation  errors,  even  though  4  bits  are  employed  in  the  adder  and 
only  3  bits  in  the  phase  shifter.  The  rms  error  for  this  case  is  0.45  ms.. 

This  computational  result  has  been  corroborated  with  an  approximate 
theoretical  analysis. 

Note  that  the  antisymmetry  property  of  t  ^  is  preserved.  Furthermore, 
it  has  been  observed  that  the  error  *  ^  present  case  with  4-bit  rcvv  and 

column  truncation  is  almost  periodic,  with  a  period  which  is  half  the  period 
given  by  the  3 -bit  phase  shifter.  The  error  is  not  exactly  periodic  at  the  half¬ 
period  intervals  because  of  the  presence  of  the  error  due  to  the  3 -bit  phase 
shifter. 

C.  Effect  of  Randomized  Quantization  Levels 

As  discussed  in  Sec,  II -D,  quantization  level  randomization  is  an 
effective  technique  for  eliminating  spatially  correlated  phase  errors  which 
give  rise  to  quantization  lobes,  i.e,,  large  angle  errors,  in  certain  directions. 
Figure  15  illustrates  the  error  obtained  in  the  absence  of  space  feed  collima- 
tion  but  with  quantization  levels  randomized.  Note  the  absence  of  any  symmetry.. 
As  expected,  the  rms  error  is  very  close  to  the  theoretically  predicted  value 
of  0.  149  ms.  The  autocorrelation  function  width  is  not  materially  different, 
however,  when  compared  to  the  autocorrelation  functions  of  the  previous 
cases.  Figure  16  show's  €^(u,  v)  in  the  same  scanning  sector  as  Figs.  8  and  11., 
Some  of  the  regular  featvxres  discernible  in  Fig.  11  are  absent  in  the  present 
case..  The  t  surface  has  a  generally  similar  character.  Figure  17  shows 
th«-  autocorrelation  function  R^^(Au,  Av).  R^(Au,  Av)  has  a  sim.ilar  shape. 

I  l)e  shape  of  tliese  surfaces  near  the  central  peak  is  very  similar  to  that 
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STEERING  ANGLE,  SINE  ALPHA  (ms)  DIFF.  STEERING 

ANGLE  (ms) 


14  Sine  Alpha  Erior  with  Low  Adder  Bit  Level  in  the  Presence 
of  Colliniation 
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STEERING  ANGLE,  SINE  ALPHA  (ms)  DIFF.  STEERING 

ANGLE  (ms) 


Fig.  17  Perspective  View  of  the  Sine  Alph.i  Autocorrelation  Function 
in  the  Presence  of  Ranrlomi/.»-cI  Quantization  Levels 
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shown  in  F’ig.  12  for  the  space  feed  collimation  case.  This  result  indicates 
that  the  spatial  statistics  of  the  angular  error  in  the  case  of  randomized 
quantization  levels  is  generally  similar  to  the  case  where  space  feed  collima¬ 
tion  is  employed.  The  width  of  the  central  peak  at  the  3  dB  point  is  still  about 
2  ms,  as  it  is  in  those  cases  employing  collimation. 

D.  Differential  Phase  Errors 

Thus  far  differenti  j.xase  errors  have  been  ignored.  Figure  18 

illustrates  the  effects  of  including  moderately  large  differential  phase  errors 

in  an  array  with  randomized  quantization  levels.  In  this  configuration  a  10°, 

one -sigma  random  differential  phase  error  is  generated  for  each  bit  with  the 

aid  of  a  gaussian  random  number  generator  and  stored  for  later  use  by  the 

computer  program.  The  errors  so  generated  are  independent  from  bit  to  bit, 

even  within  the  same  phase  shifter  element.  The  increase  in  the  rms  error 

shown  in  Fig.  18  is  representative  of  the  large  sector  average  increase  in 

rms  error  due  to  10*,  one-sigma  differential  phase  errors.  To  compare 

the  relative  magnitude  of  these  differential  phase  errors  to  the  quantization 

error  for  the  case  of  3 -bit  phase  shifters,  consider  the  total  differential 

phase  error  L  at  element  mn: 

^  ^mn 


where 


=  ith  bit  differential  phase  error,  i  =  1,  2,  3 

v.  =  ith  bit  state  function,  i  =  1,  2,  3 
1 

II  if  bit  i  is  "on" 

0  if  bit  i  is  "off" 
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STEERING  ANGLE,  SINE  ALPHA  (ms)  DIFF.  STEERING 

ANGLE  (ms) 


The  subscript  mn  has  been  omitted  from  c  .  and  V.  for  brevity.  By  appropriate 
averagi  .g  over  ihe  states  of  Yj,  and  ensemble  averaging  over  the  states  of  € 
it  is  readily  shown  that 


2  1 
“  4 


With  r,  -  (j,  =  a,  =  10°,  the  rms  differential  phase  error,  denoted  is 

*  1  ‘2  '  3  ^ 

given  by 


(r^  =  8.67*  . 

The  error  (in  the  absence  of  differential  phase  errors)  is  uniformly 
distributed  over  the  interval  ±Q/2  (=  ±22.  5°).  Therefore,  for  q  =  3  and  =  0, 

=  13.0’  . 

The  inclusion  of  the  differential  phase  error  of  8.67°  causes  the  error  o-^  to 
increase  from  13.0  to  15.6°.  This  increase  affects  cr^  and  cr^  proportionately, 
as  shown  by  Eqs.  (9)  and  (10), 

Since  the  differential  phase  errors  can  be  made  significantly  less  than 
10°  rms,  these  results  provide  good  justification  for  ignoring  differential 
phase  errors  as  an  important  source  of  angular  measurement  error.  It  is  for 
this  reason  that  differential  phase  errors  have  been  ignored  in  the  other  cases 
presented. 

A  plot  similar  to  Fig.  18  for  the  case  of  space  feed  collimation  shows 
that  the  presence  of  differential  phase  errors  destroys  the  symmetry  and 
periodicity  properties  noted  earlier.  To  the  extent  that  the  differential  phase 
errors  are  small,  these  symmetry  and  periodicity  properties  are  approxi¬ 
mately  obeyed. 
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V.  UiiJCORRELATION  TECHNIQUES  AND 
COMPUTATIONAL  RESULTS 


If  the  pulse -to -pulse  correlations  in  the  angle  measurement  error 
prove  to  be  of  consequence,  it  may  be  desirable  to  implement  one  of  several 
dithering  techniques  to  reduce  these  correlations.  The  principal  techniques 
for  accomplishing  decorrelation  c.re  (a)  phase  dithering,  (b)  frequency 
dithering,  and  (c)  beam  dithering.  These  techniques  produce  decorrelation 
of  the  angle  errors  by  forcing  random  fluctuation  of  the  element  phase 
commands  from  look  to  look. 

The  angle  correlation  is  due  to  phase  error  (4^^)  correlation  from 
look  to  look.  The  purpose  of  the  dithering  techniques,  therefore,  is  to 
generate  commands  that  vdll  tend  to  make  the  element  phase  error  uncorre¬ 
lated  from  look  to  look.  Typicall/,  the  element  phase  error  need  only 
change  by  a  phase  of  Q/4  (Q  =  2it/2*^)  to  effect  nearly  complete  decorrelation. 

A.  Phase  Dithering 

The  phase  dithering  technique  consists  of  generating  a  random  phase 
angle  (from  0  to  Zir)  and  adding  it  to  each  phase  shifter  phase  command. 

This  is  most  easily  dene  in  the  beam  steering  computer,  where  the  randomly 
generated  angle  can  be  added  to  each  row  or  column  command.  The 
resultant  change  in  the  absolut^e  phase  of  the  wave  front  leaves  the  direction 
of  the  beam  unchanged  and  is  of  no  consequence  for  many  applications.  In 
general,  the  only  situation  of  interest  in  which  the  random  phase  could  be 
undesirable  is  the  case  of  second  time  around  returns  (or  higher  order 
range  ambiguities)  from  clutter,  if  phase  coherent  clutter  rejection 
techniques  are  required. 

Figure  19  shows  the  sine  alpha  error  surface  in  the  same  scan 
sector  as  the  previous  perspective  plots.  In  this  case,  however,  phase 
dithering  has  been  added  to  the  ordinary  row  and  column  collimation 
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Collimation  and  Phase  Dithering 


configuration  whose  results  were  presented  in  Sec.  IV-B.  While  this  error 
surface  appears  to  be  much  "hashier”  than  the  surface  of  Fig.  11,  the 
really  dramatic  change  in  the  correlation  properties  can  best  be  seen  in  the 
subsequent  plot  of  R^(Au,  Av)  in  Fig.  20.  The  behavior  of  R^(Au,  Av)  is 
identical  to  R^(Au,  Av).  The  error  decorrelates  completely  in  a  differential 
angle  equal  to  the  spacing  between  locks,  in  this  case,  0.  5  ms.  Actually, 
it  is  not  necessary  to  move  the  beam  at  all  to  obtain  decorrelation  with 
phase  dithering. 

The  time -domain  autocorrelation  function  of  the  measurement  error, 
with  or  without  beam  motion,  contains  an  identically  large  central  peak  with 
complete  decorrelation  in  a  time  equivalent  to  the  inter -pulse  period,  i,  e.  , 
ihe  rate  at  which  the  phase  shifters  are  dithered.  This  technique  is  equally 
effective  in  space-fed  and  corporate-fed  arrays. 

B.  Frequency  Dithering 

The  essential  array  design  requirement  for  effective  decorrelation 
through  frequency  dithering  is  that  the  phase  shifter  phase  command  include 
an  angle  which  compensates  for  the  feed -to -phase  shifter  propagation  delay. 
Thus,  if  the  frequency  is  changed,  the  compensation  phase  changes,  result¬ 
ing  in  a  new  phase  error.  All  space-fed  planar  arrays  fall  into  this  category. 
A  corporate  fed  structure  with  varying  cable  lengths  to  each  array  element 
could  also  be  designed  to  employ  frequency  dithering. 

Frequency  dithering,  while  not  as  easy  to  implement  as  phase  dither¬ 
ing,  is  a  technique  which  will  often  be  available  for  other  reasons.  In 
particular,  frequency  agility  is  sometime  implemented  for  ECCM  or 
increased  detection  range  against  scintillating  targets.  In  these  cases, 
instrument  error  decorrelation  is  an  additional  benefit  to  consider.  Fre¬ 
quency  dithering  has  the  disadvantage  of  destroying  pulse -to -pulse  phase 
coherence  which  might  be  desirable  for  coherent  integration  and/or  clutter 
rejection. 
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Figure  21  shows  that  frequency  dithering  is  equally  effective  in 
decorrelating  the  error.  This  space-fed  configuration  was  scanned  over 
the  same  sector  as  Fig.  11.  The  total  range  of  frequency  dither  in  this 
case  is  10%  (±5%),  uniformly  distributed  across  the  band.  On  each  radar 
look,  a  new  operating  frequency  within  the  band  is  chosen  on  a  random 
basis.  It  is  easy  to  show  that  a  5%  change  in  frequency  causes  most  of  the 
phase  shifters  to  change  by  at  least  Q/4. 

C.  Beam  Dithering 

The  angle  errors,  as  shown  in  Sec.  IV-B,  indicate  that  correlations 
for  the  particular  array  studied  never  exceeded  more  than  a  few  millisines. 

A  natural  decorrelation  technique,  therefore,  is  to  steer  the  beam  slightly 
more  than  necessary,  in  a  random  way,  so  as  to  enhance  the  natural 
decorrelation  that  already  exists.  The  resulting  degrade tion  due  to  off-axis 
tracking  is  small,  since  the  total  excursion  due  to  dithering  can  be  chosen 
to  bo  much  less  than  a  bcamwidth.  The  implementation  of  this  technique 
is  also  quite  simple.  It  is  entirely  contained  in  the  beam  steering  computer. 
A  random  number  generator  selects  two  angles  within  specified  limits. 

These  are  then  simply  added  to  the  beam  steering  computer  input  steering 
commands. 

In  order  to  generate  an  error  surface  which  has  the  correlation 
properties  shown  in  Fig.  22,  a  random  angle  generator  was  used  which 
produced  variates  with  a  uniform  distribution  from  -5  ms  to  +5  ms,  or 
about  20%  of  the  beamwidth.  The  beam  is  dithered  to  a  new  random  angle 
(both  u  and  v)  on  each  radar  look.  The  correlation  function  is  seen  to  drop 
very  rapidly  to  a  small  pedestal  and  then  becomes  very  small.  This  tech¬ 
nique  n.ight  be  improved  by  changing  the  distribution  of  the  randomly  varying 
angle  commands,  although  this  possibility  has  not  been  investigated,  A 
disadvantage  of  beam  dithering  is  that  the  total  beam  excursion  is  limited 
somewhat  by  off-axis  degradation.  In  situations  where  MTI  is  employed  to 
reject  clutter,  beam  dithering  may  seriously  degrade  MTI  performance.  It 
is  concluded  that  beam  dithering  is  a  good  decorrelation  technique  but  is  less 
attractive  than  frequency  dithering. 
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Fig.  22  Perspective  View  of  the  Sine  Alpha  Autocorrelation  Function 
with  Collimation  and  Beam  Dithering 
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VI.  SUMMARY 


This  paper  describes  an  analysis  of  the  random  instrumentation  angle 
measurement  errors  of  phased  array  radars.  A  number  of  calculations  of 
angle  measurement  errors  were  performed  for  a  particular,  though  typical, 
array  configuration.  These  calculations  were  aided  by  a  comprehensive 
computer  program,  which  was  de /eloped  in  sufficiently  general  form  to 
permit  its  application  to  a  wide  class  of  phased  array  radars.  The  program 
simulates  finite  computation  accuracy  in  the  beam  steering  computer,  trun¬ 
cation  of  steering  commands  at  the  row  and  column  level,  truncation  of 
steering  commands  at  the  phase  shifters,  and  the  effects  ol  differential  phase 
errors.  The  example  array  was  examined  in  various  modes  of  operation, 
with  and  without  quantization  level  randomization  and  with  and  without  space - 
feed  collimation.  The  effects  of  row  and  column  adder  computation  accuracy 
were  examined.  The  statistical  behavior  of  the  angle  measurement  error 
with  and  without  several  dithering  techniques  was  determined.  For  the  typical 
array  configuration  considered,  the  following  conclusions  can  be  drawn; 

a.  A  beam  split  ratio  (neglecting  bias  errors)  in  excess  of  300:1 
can  be  achieved  with  the  use  of  only  3-bit  phase  shifters. 

b.  Truncation  of  row  and  column  steering  commands  was  found  to 
dominate  the  angle  error  in  certain  array  configurations.  In 
these  cases  the  beam  split  ratio  may  be  reduced  below  100:1 
and  the  angle  error  is  nearly  periodic  at  a  submultiple  of  the 
array  period. 

c.  Differential  phase  errors  as  large  as  lO”,  one- sigma  per 
bit  were  found  to  have  only  a  minor  effect  on  rms  angular 
measurement  error.  The  rms  error  increased  only  about 
20%.  A  bias  error  on  the  order  of  0.  1  m.s  was  also  observed. 

d.  A  steering  angle  of  approximately  2  ms  in  any  direction  will 
cause  a  50%  drop  in  the  autocorrelation  of  the  angular  errors. 
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PrecBding  p8|8  Wank 


The  horizontal  cross  section  of  the  central  peak  of  the 
two-dimensional  autocorrelation  function  is  roughly  circular, 
i.e.  ,  there  are  not  strong  "preferred"  correlation  directions. 

e.  The  phase  dithering  technique  was  shown  to  be  completely 
effective  in  annihilating  all  angle  measurement  error  correla¬ 
tions.  Results  were  presented  for  a  typical  sector  scan,  but 
beam  steering  is  not  necessary.  This  technique  is  very  easy 
to  implement. 

f.  Frequency  dithering  was  also  shown  to  be  an  effective  decorrela 
tion  technique.  Beam  dithering  was  found  to  be  somewhat  less 
effective  than  frequency  or  phase  dithering.  Both  frequency  and 
beam  dithering  involve  compromises  which  may  make  them 
unacceptable  for  certain  applications. 


54 


APPENDIX  A 

RMS  ERROR  FOR  VARIOUS  APERTURE  WEIGHTING  FUNCTIONS 


The  Eqs.  (9)  and  (10)  show  th«ii  me  rms  angle  measurement  error  is 
dependent  on  the  diffcxence  pattern  amplitude  weighting  functions  6^  and  6^. 

It  is  the  purpose  of  this  appendix  to  show  that  the  rms  error  is  only  weakly 
dependent  on  the  actual  shape  of  the  difference  pattern.  It  is  intuitively 
obvious  that  those  weighting  functions  that  tend  to  place  more  weight  on  the 
phase  shifters  at  the  edges  of  the  array  will  tend  to  produce  lower  rms  angle 
errors.  On  the  other  hand,  such  functions  tend  to  produce  larger  sidelobes. 

To  illustrate  the  effects  of  various  difference  pattern  weighting  functions  on 
the  rms  error,  Eqs.  (9)  and  (10)  were  numerically  evaluated  on  a  small 
programmable  desk  calculator  (HP910tA).  The  results  are  compared  in 
Table  A-1. 

Only  the  x-difference  pattern  weighting  rmction  is  illustrated  in 
Table  A-1.  Note  that  this  is  assumed  to  be  independent  of  y.  The  y-differencc 
pattern  weighting  function  was  assumed  to  be  the  same  as  the  x-difference 
pattern  with  the  variable  y  replacing  x.  To  summarize, 

6x(*»y)  = 

6y(x,y)  =  6y(y)  = 

Other  applicable  parameters  are  as  defined  in  Sec.  IV.  The  results  in 
Table  A-1  show  that  the  rms  error  spread  is  likely  to  be  less  than  25%  of  a 
representative  mean  value  for  all  those  weighting  functions  of  practical 
interest. 
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Table  A-1  RMS  Angle  Measurement  Er-.or  for  Various 
Amplitude  Weighting  Functions 
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